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The volume common to two cylinders of radii r l <r 2 with axes intersecting at angle 
|8 is found to be rlv(k)/sm /3, where k=rjr 2 and v(k) may be evaluated (1) as the hypergeo- 
metric series 



2^F I (-U;2;*°)=4*$(*)( n i 1 )**. 



(2) as the combination of complete elliptic integrals (8/3)[(l-\-k 2 )E(k)- 

(3) as the cumulative integral 



■(l-k^Kik)] or 



j; 



8 kE(k)dk. 



A table of v(k) to 8 decimals over the range 0< ft (0.01) <1.00, including 5 2 m modified second 
central differences, is presented. This volume integral was useful in interpreting a gamma- 
ray albedo experiment involving a collimated source and a collimated detector, and may also 
be applicable to crossed-beam experiments. Two series useful for k close to unity are pro- 
vided, one of which involves differencing against the series 



16/3=4. $(i)( B -i> 



1. Introduction 

In crossed-beam experiments [1] 1 using the high- 
intensity accelerators now becoming available, the 
"geometrical target/' or volume common to the two 
colliding beams, is a useful parameter for interpreting 
the measured data. An evaluation of this volume in 
terms of an infinite series was recently exhumed for 
possible application to an x-ray free-air ionization 
chamber having a cylindrical sensitive volume inter- 
sected by a pencil of x rays [2]. This evaluation 
had been used in the analysis of a gamma-ray beam 
back-scattering experiment [3] for making a theoreti- 
cal estimate of the single-scattered component of 
the radiation "seen" by a collimated detector. 

Evaluations of the volume common to two circular 
cylinders of unequal radii with axes intersecting 
at right angles [4, 5, 6], and of equal radii with axes 
intersecting at an arbitrary angle [7], have frequently 
been offered as calculus textbook exercises. How- 
ever, a combined treatment does not seem to appear 
in the technical literature in a form convenient for 
easy application to practical problems. The follow- 
ing results provide formulas, a table, and a graph 
for such applications. 

The series used in [3] is here corrected, expressed 
in terms of binomial coefficients, and identified as a 
hypergeometric series. For nearly equal cylinder 
radii, convergence can be accelerated by use of the 
difference-series technique [8]. An alternative series 



i Figures in brackets indicate the literature references at the end of this paper. 



for this region, derived from the right-angle elliptic- 
integral solution [4, 5, 6], is more complicated but 
also more rapidly convergent. 

2. Volume Integral 

The integral for the common volume of two 
cylinders of radii ri<r 2 with axes intersecting at 
angle /3 (see fig. 1) is found as follows. The cross 
section parallel to the cylinder axes, at a distance 
x from them, is a parallelogram of height 2{r%— x 2 ) 1/2 
and base 2{r\— x 2 ) I/2 /sin fi. Hence the volume inte- 
gral is 

2{rl-x*y» - Kn . X> dx (la) 
-rj sin p 

=~ f ri (ri-x*y< 2 ■ (rl-x 2 y'*dx. (lb) 

sin p Jo 

3. Common Volume When r x =r 2 

For equal cylinder radii r 1 =r 2 =r, the integral in 
(lb) reduces to the familiar result [7] 



sinp jo 



16^ 



3sin/3 



(2) 
(3) 
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Figure 1. Three-view sketch of the common volume of cylinders 
with radii r t and r 2 axially intersecting at angle ft. 

The area of the shaded parallelogram parallel to the plane of the axes in the lower 
left view comprises the integrand in eq (la) and is integrated over the range 
—r\<x<n shown in the other two views. 

4. Common Volume When r l <r 2 

4.1. Series Solution 

The factor (r%— x 2 ) 1 ' 2 in the integral in (lb) may 
be expanded as a power series [9, p. 2, eq 5.3] 
in x/r 2 , since x<r x <r 2 . The volume integral then 
becomes 

which may be integrated term by term. 
The resulting series solution is 

_3 / 1 v AiY-5 ( l ' s Y ( r A* 

S\2-4j \rj 8\2-4-6/ W 

1 ( 1 - 3 - 5 V friY - \ 

10 \2 ■ 4 • 6 • 8/ W " ) 

= Ss(»)C-i)w ■ (5) 



This same result may be obtained by casting the- 
integral in (lb) in the form 

F(r 1) r 2 ,j8) = (^/sin/3)4(^y 

x J l r l ' 2 (i-ty 2 (i-^| tj' 2 dt, (6> 

where t=x 2 frl, which is now recognizable as the 
integral representation of the hypergeometric series 
[10] 

V{r u r 2 , « = (rf/sin /3)2x (£)>, (-1 \, 2; J) (7) 

identifiable with (5). 

Since (5) is somewhat slowly convergent when 
r x ^r 2 , under some circumstances it may be advan- 
tageous to difference this series against a I/71-- 
series (16) discussed in the appendix, giving 

(8) 



V(r u r h p) = 



4wrl 
sin/ 



The convergence rate of the series-term in (8) is 
not improved over that of (5). However, for 
r l ^r 2 this sum is small compared to the constant 
term 4/St, identifiable with the equal-radii solution 
(3), hence resulting in higher precision of V(r lf 
r 2 , /3) for the same number of terms. An alterna- 
tive series solution for this region is given at the 
end of the following section. 

4.2. Elliptic Integral Solution 

An alternative solution of the integral in (lb) may 
be obtained as a combination of complete elliptic 
integrals [11] of the first and second kinds, K(k) 
and E(k). Applying formulas (219.11) and (361.03) 
from Byrd and Friedman [12] the result is found to 
be 2 



V(r u r 2t P)- 



Srl 



3sinj3 



where 



[{l+k 2 )E{k)-{l-V)K{k)\ (9) 



k=ri[r 2 . 



Except for the angle factor 1/sin /3 this result is 
the standard textbook solution [5, 6] for cylinders 
intersecting at right angles. Also, the formulation 
in (9) is related to the indefinite integral [12, eq 
(611.01)] 



x 



* kE{k)dk=\[(\- 



-¥)E{k)-{l-k 2 )K{k)}. (10) 



2 This integral is part of the "G factor" used for interpreting gas scattering ex- 
periments in which a circular-aperture detector views a gas target transversed 
by a cylindrical beam. In this context this elliptic integral solution has been 
given by E. A. Silverstein, Nucl. Instr. and Meth. 4, 53 (1959) and by D. F. 
Herring and K. W. Jones, Nucl. Instr. and Meth. 30, 88 (1964). 
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A partial check on (9) is obtained by expanding 
E(k) and K{k) as power series in k according to 
reference [9], eqs (773.1) and (774.1). Combining 
like powers of k and substituting back r x \r 2 for k the 
resulting series is identical with (5). An additional 
check is provided by the identity of eqs (9) and (3) 
in the limit as &->l. Also, eq (9) can be obtained 
from (7) by use of the tenth Gauss recursion formula 
on page 9 of reference [10]. 

For k close to unity, a series which converges more 
rapidly than (8) may now be derived by substituting 
in (9) the series in equations (773.3) and (774.3) in 
reference [9] for K{k) and E{k). The first few terms 
of this series are 

F(r„*,/3)=W/sin p) -y { 1-f k' 2 

2MV V V 2 2 -4 2 -6V *' 12/ 

l 2 -3 2 -5 A £_43\,, 8 
d '2 2 -4 2 -6 2 -8\ k' 40/ 

. l 2 -,3 2 - 5 2 -7 /, 4 967\,,io \ ,„. 

- 4 >-4^6»-8M0 V ln F~840r -•••/ (U) 

where k' 2 =l— k 2 = l — (ri/r 2 ) 2 . Using only the terms 
given in (11) the sum for ri/r 2 =0.90, without the 
factor (rg/sin/S), gives 4.49991482 as compared 
with the exact value of 4.49991288 . . ., and the 
convergence improves as r x \ri goes toward unity. 
The series in (11) may be obtained in general form, 
if desired, by use of the appropriate transformation 
[13] on the hypergeometric series given in (7). 

5. Numerical Results 
In table 1, the dimensionless factor 

sin /3 



v(k)-- 



r\ 



V(r lt r„ 0), 



-*£®Cii)*"' 

=8 C kE(k)dk, 



(12a) 

(12b) 
(12c) 
(12d) 



=U(l+k 2 )E(k)-(l-k 2 )K(/c)], (12e) 

where k=rjr 2 is tabulated to 8 decimal places for 
0<fc(0.01)<1.00, computed using (12e) and K(k) 
and E(k) from [11]. Hence, for many practical 
applications, the common volume of two cylinders 
with radii r^^ and axes intersecting at angle /3 
may be computed as 



v^n-Kp 9 ® 



(13) 



in which values of v{k) are interpolated from table 1. 

Table 1. Values of v(k), defined in eq (12a-e), over the range 
0<k(0.01) <1.00, valid to the 8D given 

Modified second central differences 8 m are provided for interpolation using 
auxiliary tables. 



ft 


v(k) 


*i 


k 


v(k) 


4 


0.00 


().()()()() 


0000 


+12 


5665 












.01 


0.0006 


2831 


12 


5656 


0.51 


1.5792 


3625 


+9 


8789 


. 02 


0.0026 


1315 


12 


5627 


.52 


1 . 6394 


2033 


9 


7607 


.03 


0.0066 


5423 


12 


5580 


. 53 


1.7005 


8043 


9 


6391 


.04 


0.0100 


5109 


12 


5514 


.54 


1.7627 


0439 


9 


5142 


.05 


0.0157 


0305 


12 


5429 


.55 


1.8257 


7969 


9 


3854 


.06 


0. 0226 


0928 


12 


5325 


.56 


1.8897 


9346 


9 


2530 


.07 


0.0307 


6874 


12 


5203 


.57 


1 . 9547 


3247 


9 


1167 


.08 


0.0401 


8019 


12 


5061 


.58 


2. 0205 


8308 


8 


9766 


.09 


0.0508 


4222 


12 


4900 


.59 


2. 0873 


3128 


8 


8323 


.10 


0. 0627 


5321 


12 


4720 


.60 


2. 1549 


6262 


8 


6837 


.11 


0. 0759 


1138 


12 


4520 


.61 


2. 2234 


6226 


8 


5308 


.12 


0. 0903 


1471 


12 


4302 


.62 


2. 2928 


1490 


8 


3734 


.13 


0. 1059 


6104 


12 


4065 


.63 


2. 3630 


0478 


8 


2112 


.14 


0. 1228 


4797 


12 


3807 


.64 


2. 4340 


1570 


8 


0441 


.15 


0. 1409 


7294 


12 


3530 


.65 


2. 5058 


3093 


7 


8719 


.16 


0. 1603 


3316 


12 


3234 


.66 


2. 5784 


3325 


7 


6943 


.17 


0. 1809 


2569 


12 


2917 


.67 


2. 6518 


0400 


7 


5113 


.18 


0. 2027 


4735 


12 


2581 


.68 


2. 7259 


2757 


7 


3222 


. 1!) 


0. 2257 


9478 


12 


2224 


.69 


2. 8007 


8235 


7 


1272 


.20 


0. 2500 


6441 


12 


1848 


.70 


2. 8763 


4974 


6 


9258 


.21 


0. 2755 


5248 


12 


1451 


.71 


2. 9526 


0957 


6 


7174 


. 22 


0. 3022 


5502 


12 


1034 


.72 


3. 0295 


4102 


6 


5020 


'.23 


0.3301 


678(5 


12 


0595 


.73 


3. 1071 


2253 


6 


2791 


.24 


0. 3592 


8660 


12 


0136 


.74 


3. 1853 


3179 


6 


0481 


.25 


0. 3896 


0667 


11 


9656 


.75 


3. 2641 


4571 


5 


8086 


.26 


0.4211 


2326 


11 


9154 


.76 


3. 3435 


4031 


5 


5600 


.27 


0.4538 


3135 


11 


8632 


.77 


3. 4234 


9073 


5 


3016 


.28 


0. 4877 


2572 


11 


8087 


.78 


3. 5039 


7111 


5 


0328 


.29 


0. 5228 


0002 


11 


7520 


.79 


3. 5849 


5457 


4 


7527 


.30 


0. 5590 


5129 


11 


6932 


.80 


3. 6664 


1308 


4 


4605 


.31 


0. 5964 


7093 


11 


6320 


.81 


3. 7483 


1740 


4 


1550 


.32 


0. 6350 


5373 


11 


5686 


.82 


3. 8306 


3696 


3 


8350 


.33 


0. 6747 


9336 


11 


5020 


.83 


3. 9133 


3972 


3 


4992 


.34 


0. 7156 


8323 


11 


4348 


.84 


3. 9963 


9206 


3 


1456 


.35 


0. 7577 


1654 


11 


3644 


.85 


4. 0797 


5861 


2 


7725 


.36 


0. 8008 


8626 


11 


2915 


.86 


4. 1634 


0199 


2 


3774 


.37 


0. 8451 


8508 


11 


2163 


.87 


4. 2472 


8265 


1 


9574 


.38 


0. 8906 


0549 


11 


1385 


.88 


4.3313 


5853 


1 


5089 


.39 


0. 9371 


3970 


11 


0582 


.89 


4.4155 


8469 


1 


0276 


.40 


0. 9847 


7969 


10 


9753 


.90 


4. 4999 


1288 


+ 


5075 


.41 


1.0335 


1716 


10 


8899 


.91 


4. 5842 


00! Hi 


_ 


586 


.42 


I.0N33 


4356 


10 


8017 


. 92 


4. 6686 


6213 




6811 


.43 


1. 1342 


5009 


10 


7108 


.93 


4. 7529 


6390 


1 


3730 


.44 


1. 1862 


2764 


10 


6171 


.94 


4. 8371 


2665 


2 


1554 


.45 


1.2392 


6686 


10 


5207 


.95 


4. 9210 


7158 


3 


0581 


.46 


1. 2933 


5809 


10 


4213 


. 96 


5. 0047 


0740 


4 


1308 


.47 


1. 3484 


9141 


10 


3191 


.97 


5. 0879 


2495 


5 


4606 


.48 


1. 4046 


5657 


10 


2137 


.98 


5. 1705 


8678 


-7 


1824 


.49 


1.4618 


4304 


10 


1053 


.99 


5. 2525 


0293 







.50 


1. 5200 


3999 


+9 


9937 


1.00 


5. 3333 


3333 








Modified second central differences b 2 m are included 
for interpolation by Everett's formula 



v(k i +pAk) = (l-p)v(k i )+pv(k i+l )+E 2 (p)5 2 m , i 

+F 2 (p)8l, i+1 + . . 



(14) 



where p is the interpolation fraction of the interval of 
tabulation Ak, and 8 2 mti) 5^ )Z+ i are the modified second 
differences at the tabular points i and i+1 and were 
evaluated from the second and fourth differences 
according to 

dl=d 2 — 0.1845 4 . 
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Everett's coefficients E 2 (p) and F 2 (p) are available 
in standard tables [14, 15] and are identical with 
Lagrangian interpolation coefficients A A ^ and At 
[16, table 25.1]. 

The general behavior of v(k) is shown by the curve 
in figure 2. Values of v(k) can be taken directly 
from this curve for use in rough calculations where 
only two- or three-figure accuracy is required. 



v(k) 
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Figure 2. Graph of v(k) showing the general behavior of the 
function and suitable for rough calculations. 



6. Appendix. Two Series for 1/tt 

The series evaluation (5) for V(r 1 ,r 2 ,P) contains a 
factor of x and the formula (3) for V{r,$) does not. 
Thus, for r 1 =r 2 =r, the right-hand side of (5) can 
be equated to the right-hand side of (3) to form 



16r 3 



3sm/3 



=sSsC)U0 (i5) 



from which 



3tt £=[\nj\n-lj 2 4\2/ 6V2-4/ 

5/ 1-3 V 

8 \2 • 4 • 6/ 



(16) 



This series can now be used to form the difference- 
series in eq (8). 



An additional l/V-series, which also does not appear 
in standard compilations of series [17], can be ob- 
tained by combining (16) with a series discussed by 
Bromwich [18], [17, eq 274] 



1 » / 1 \2 

-=yw 2 ) 

r n =i\n— 1/ 



+ 



to form 






Using the addition theorem 

GHZH'i 1 ) 



(17) 
(18) 

(19) 



the result is 



3tt n=i\n—\J\n/ 



-1 5+i z ' 1 I 1 3-1-1 

*2 i "2'2.4" i "2.4*2-4-6 

1-3 3-1-1- 3 



2-4-6 2-4-6-8 
1-3-5 3-1-1- 3-5 



2-4-6-8 2-4-6-8. 10 



(20) 
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the problem, and P. Lamperti and L. W. B. Jolley 
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section and 1/tt series, respectively. The author is 
also indebted to F. W. J. Olver, I. A. Stegun, L. F. 
Epstein, and A. Fletcher for suggestions and com- 
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computing table 1 . 
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